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Проведено аналіз сучасних підходів до пошуку елементів на зображеннях. Розроблено алгоритм пошуку 
мікрооб'єктів на гістологічних 1 цитологічних зображеннях з використанням бази даних. Розроблено ярусно- 
паралельну форму розпаралелення процесу пошуку мікрооб'єктів за шаблоном. Здійснено програмну реаліза- 
цію пошуку мікрооб'єктів за шаблоном. Отримані результати демонструють, що час роботи програмного мо- 
дуля із застосуванням процесу розпаралелення пришвидшує процес опрацювання в середньому на 2096 для 
цитологічних зображень. 
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Вступ 

На даний час у медицині для діаг- 
ностування ракових та передракових ста- 
нів молочної залози широко використо- 
вують гістологічний і цитологічний аналіз 
ПІ). Із розвитком апаратних та програмних 
засобів, із застосування засобів штучного 
інтелекту та комп'ютерного зору, значно 
збільшилась кількість програмних засобів 
для аналізу зображень. Використання су- 
часних ШП-технологій значно підвищує 
продуктивність та якість роботи лікарів- 
діагностів. 

Форма та взаєморозміщення ядер 
клітин на цитологічних зображеннях може 
вказувати на наявність певних патологій, 
які потребують детальнішого досліджен- 
ня. Однак, зображення, отримані в резуль- 
таті мікроскопічного дослідження клітин, 
характеризуються високим рівнем різного 
роду шумів, що ускладнює процес автома- 
тичної обробки та виділення ядер клітин. 


До того ж, ядра клітин у більшості випад- 
ків не є чітко виділеними, а навпаки - роз- 
митими та такими, що дотикаються одне 
до одного. 

Однією з технологій пошуку елемен- 
тів на зображенні є пошук елементів за на- 
перед заданим шаблоном. Пошук елемен- 
тів за шаблоном - це техніка цифрової об- 
робки зображень для пошуку невеликих 
частин зображення, які відповідають зада- 
ному зразку. У більшості випадків лікар 
обробляє стек цитологічних зображень, 
отриманих з камери мікроскопа. Кожне 
зображення характеризується з високою 
роздільною здатністю, тому пошук навіть 
декількох наперед заданих шаблонів на 
вхідних зображеннях може значно заван- 
тажити робочу станцію лікаря-діагноста. 

Для пришвидшення процесу пошуку 
елементів за шаблоном у даній статті роз- 
роблено основні кроки розпаралелення да- 
ного процесу. 
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Постановка проблеми 

Отже, проблемою, що існує, є роз- 
роблення та впровадження засобів пошуку 
елементів на зображенні за наперед зада- 
ним шаблоном та пришвидшення процесу 
їх обчислення, формування бази шаблонів 
патологічних станів. Це дозволить значно 
підвищити якість досліджень та спростити 
процес постановки діагнозу. 

Аналіз публікацій 

Науковці значну увагу приділяють 
автоматизації процесу діагностування 
передракових та ракових станів. Оскільки 
одним із основних способів діагносту- 
вання онкологічних захворювань є аналіз 
зображень зрізів тканин, то значна увага 
приділяється засобам розпізнавання та 
класифікації зображень. Наприклад, авто- 
ри 12,3| пропонують алгоритм виділення 
ядер клітин та мікроскопічних зображень 
на основі їх форми та забарвлення. 

У роботі (4) наведено метод пошуку 
об'єктів на зображенні, базуючись на 
афінних перетвореннях та максимальній 
точці кореляції для подальшого макси- 
мального накладання. 

У статтях |5,6| наведено аналіз засо- 
бів для розпаралелення процесу виділення 
контурів ядер клітин на широкоформат- 
них зображеннях за допомогою графічних 
процесорів. Автори у статті |7| наводять 
порівняльний аналіз процесу опрацю- 
вання біомедичних зображень із та без ви- 
користання розпаралелення на графічних 
процесорах, що доводить актуальність за- 


Модуль завантаження стеку 
зображень 


Фільтрація (медіанний, 
гаусовий) фільтри 


Модуль аналізу зображень 
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стосування різних підходів до розпарале- 
лення обчислень для пришвидшення про- 
цесу діагностування. 

Важливим етапом у процесі пошуку 
об'єктів за шаблоном є підготовка зо- 
браження - попереднє оброблення, філь- 
трація, гістограмне вирівнювання. Автори 
Реїг 5шгако5 та Міап /аго5 у роботі |8,9| 
наводять різні підходи та результати роз- 
паралелення процесу фільтрації біомедич- 
них зображень. 

У статтях |10,11| наведено порів- 
няльний аналіз існуючих  програмно- 
апаратних засобів для автоматизації про- 
цесу діагностування гістологічних та ци- 
тологічних зображень. Наявність алгорит- 
мів розпаралелення є значною перевагою. 

Мета дослідження 

Отже, метою роботи є розробка уза- 
гальненої структури модуля пошуку еле- 
ментів на цитологічному зображенні за 
шаблоном та його розпаралелення для 
системи «НІАМ5» (12| для підвищення 
швидкості та точності діагностування ра- 
кових захворювань. 

Узагальнена структура модуля 
аналізу зображень 

Модуль аналізу зображень є основ- 
ною частиною програмної системи, який 
здійснює такі функції: завантаження стеку 
вхідних зображень у пам'ять, попередньої 
обробки, фільтрації, завантаження напе- 
ред заданих шаблонів, пошук об'єктів за 
ними і вивід кінцевих результатів. Струк- 
тура модуля наведена на рисунку |. 


Модуль зберігання 
шаблонів 


Модуль пошуку об'єктів на 
вхідних зображеннях 


Рис. 1. Структура модуля аналізу зображень 
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Модуль завантаження у пам'ять 
отримує файли вхідних зображень, на 
яких буде проводитися пошук об'єктів, та 
завантажує їх у базу даних програми. За- 
вантаження зображень проводиться як із 
внутрішнього сховища комп'ютера, так і з 
зовнішніх носіїв чи мережі Інтернет. 

Після завантаження вхідних зобра- 
жень відбувається їхня попередня оброб- 
ка, а саме: фільтрація медіанним або гау- 
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Гістограмне 
вирівнювання 


Фільтрація 


Завантаження 
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Завантаження 
шаблону 


Зображення 


Гістограмне 
вирівнювання 


Фільтрація 


совим фільтром. Небажаний ефект від за- 
вад усувається фільтрацією. 

Розпаралелення процесу пошуку 
об'єктів на зображенні 

Пошук за шаблоном на декількох зо- 
браженнях є незалежною задачею, тому 
доцільно застосовувати процес розпарале- 
лення. Ярусно-паралельну форму пред- 
ставлення розпаралелення процесу пошу- 
ку елементів за шаблоном наведено на 


рисунку 2. 


Зображення 
ЯМ 


Гістограмне 
вирівнювання 


Накладання Накладання Накладання 
шаблону на шаблону на шаблону на 
зображення зображення зображення 
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результату результату результату 


Рис. 2. Ярусно-паралельна форма представлення 
розпаралелення процесу пошуку елементів за шаблоном 
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На першому етапі відбувається за- 
вантаження шуканого шаблону в пам'ять 
комп'ютера. На наступному - проходить 
розгалуження процесу обробки зобра- 
жень, оскільки дані етапи є незалежними 
і, відповідно, їх можна розподілити на 
окремі потоки. У результаті роботи на ви- 
ході отримується нове зображення і, якщо 
на ньому присутні подібні ядра клітин, то 
вони будуть відповідно виділені. Надалі 
лікар-діагност матиме можливість деталь- 
ніше звернути увагу на ті цитологічні 30- 
браження, на яких знайдено відповідні 
області. 

Алгоритм пошуку мікрооб'єктів 
за шаблоном 

Алгоритм пошуку мікрооб'єктів за 
шаблоном складається з таких кроків: 

1) Виділення контурів ядер клітин За ал- 
горитмом, розробленим у статті | 131. 

2) Ідентифікація трьох точок на контурах 
ядер клітин і обчислення афінних кое- 
фіцієнтів 114). 
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Коефіцієнти афінних перетворень, 
що описують зміну форми і розміщення 
початкового ядра, отримуються із системи 
рівнянь: 


дл у 31 Та Ж 
х, у; 1||5|«|Х, 
х» ЗУ; Ще х 


д 
Д 
інь 
чаї 


, 


де (хі, у) й (Ха, Уз) , (Хз» Уз) - точки, що 
належать першому ядру; 


(ам), (з), (хз) точки, що нале- 
жать другому ядру; 

а, Б, с, 4, е, Ї - коефіцієнти афінного пере- 
творення, що відображають перше ядро 1 


друге. 
Невідомі коефіцієнти визначаються так: 


Рани у (8 25) ж уз 2 о) уз (2 - Ху) 


(з зх.) ж У (ху 2 0) ж Уз(Х» зх) 


Б 


, 


- хі б охо (б о ху) (0, о Хо) 


(б з Хо) з (ху з Ха) У» зх) , 


рр хі(уз35 7 Узхз Хо (Уа; з ДЗК (Хо з Уз ху) 


(3 з 22) У (ху з ха) Ж У (3, з Ху) 


да у (з 7 У») ж уз (У) з Уа) Ж Уз (У» - У) 


(з з ха) Ж У (о п) У» зх) б 


я хі (95 з Узд охо (Уа з У) ха (У; - У») 


(з ох.) ж У (ху 25) Ж Уз (0, зх) 


, 


Зо хі(узУз -узУг) ях» (уз - ув) (У У» - угУ) 
у (а з 5) У (ху з Ха) Ж Уз(Х» з х,) 


3) Пересування шаблону по усіх можли- 
вих позиціях та розрахунок для кожної 
позиції «подібності» між шаблоном 1 
областю на зображенні. 

4) Обчислення значення найбільшої подіб- 
ності здійснюється згідно з формулою: 
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де Її - значення пікселя шаблонного 
зображення, 

Х - середнє значення пікселів у шаб- 
лонному зображенні, 

У | значення пікселя вхідного 
зображення, 

У || середнє означення пікселів 


вхідного зображення, 
М -- загальна кількість пікселів. 


Комп'ютерні експерименти 

Для тестування роботи розробленого 
програмного модуля було оброблено 
цитологічні зображення із бази даних 
«ВРСІ2100» (15). Вибірку зображень поді- 
лено на такі класи: (1) цито-фібро-кістозна 
мастопатія; (2) цито-кістозна мастопатія; 
(3) цито-непроліферативна фібромастопа- 
тія; (4) цито-непроліферативна мастопатія. 

Розроблений модуль є частиною 
системи | автоматизованої мікроскопії 
«НІАМ5». Результат роботи програми на- 
ведено на рисунку 3. 


Шлях до зображення 
СЛьоптогюм модії вер 
С Льовтуленоду 03 яр 
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С Льомтолемю 007 ері 13 05.2017.14.15 161 6о 
С Льомтдлем о (08 вер 

СЛьютлечних 0 вето? 05.2017. 16.05.0019 
СЛьритолемниор 008 ретрігВ 04.2017.21.53.27)9мр 


С Лькютоа ето 009 вто 


" Площа Периметр 
з 150555 447.20309817790985 
55 «вто 255 634991966847 
5 тма65 «оо 422024011612 
-? 161330 чазотутуєтмизом 
ю тотаао 131 44573859400085 
в аа 5 309 74725М1263174 
Го 1333 ча музою 
є дуб 5 мтлтовзазууттта 
т 6095 тт0.08029У920372 
16435 196.46803У9363556 
те 518175 1560.9129003286262 


Рис. 3. Графічний інтерфейс вікна системи автоматизованої мікроскопії «НІАМ5» 


На рисунку 4 наведено приклад вхід- 
ного зображення та шаблонного зобра- 


ження. Результат пошуку об'єкта наве- 
дено на рисунку 4(в). 
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а) вхідне зображення 6) шаблон 


в) результат 


Рис. 4. Пошук об'єкта на зображенні за шаблоном 


Результати роботи програмного мо- зображеннях розміром 3000 х 3000 піксе- 
дуля для пошуку об'єкта на 20-и вхідних лів наведено на рисунку 5. 
60000 56696 
45352 

45000 о рр 
о 
5 
сзобо0 
(аз) 
г 

15000 

0 
Без розпаралелення З розпаралеленням 


Рис. 5. Час пошуку об'єктів за шаблоном для стеку із 20-и зображень 
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Як видно з рисунку 5, час роботи 
програмного модуля з пошуку об'єктів за 
шаблоном на 20-ьох широкоформатних 
зображеннях приблизно на 2095 менший, 
ніж без використання розпаралелення. Бу- 
ло проведено 5 спроб для кожного випад- 
ку 1, так як час роботи програм був різний, 
на графіку відображено середній час для 
кожного випадку. Даний результат під- 
тверджує актуальність задачі розпарале- 
лення процесів опрацювання зображень. 

Висновки 

У запропонованій роботі розроблено 
алгоритм пошуку мікрооб'єктів за шабло- 
ном та ярусно-паралельну форму процесу 
пошуку мікрооб'єктів на зображенні за 
шаблоном. Отримані результати показу- 
ють, що ефективність застосування проце- 
су розпаралелення складає близько 2092 
для цитологічних зображень із бази даних 
«ВРСІ2100». 
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О.М. Вегел7зКу, О.У. Ріїзип, 
Т.М. Роїупиук, Хи.М. Ваїко 

Рага!еіїпо ої Ше ргосез55 ої 5еагсрпіпе 
обіесі5 оп суїоіовіса| ітаєеє Бу Ше 
тетріате 

Тре и5е ої піодегп ГТ іесплпоїогіе5 іп 
тедісіпе, уурісп 5ієпійсапіу іпсгеа5е Ше 
еййсіепсу апа диайсу ої (Де ФШаєпобіїсіапзя" 
ууогК, 15 апаїугед. Діавповіпе раФоїовісаї 
ргосе85е5 іп опсоіоєу 15 Бба5ед оп Ше апа- 
Туз15 ої біотедіса! (суоіовіса! апа Бізіоіо- 
дісаї) ітаєбе5. СепегаПу, сукоіобіса! апа 
різсоіовіса! ітаєє5 аге срагасіегілед Бу 
рієр сотріехіїу ргосе55іпє. ТБе піаїп 4і8- 
адуапіаєез аге Ше рге5епсе ої аіНегепі 
Кіпд5 ої пої5е апад рієп деппіцоп. Су(оЇобі- 
са! іта?єе5 аге апаЇугеа. ПІ 15 5рохуп Фаї Фе 
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5раре апа іп(егрозійоп ої сей писіеі іп 
сушіовіса! ітаєе5 тау іпдісаїе Ше рге- 
5епсе ої сегіаїп раШФоїоб1іе5. Тре іесппоїову 
ої ітаєе еіетепіз райегп 5еагсП 15 п5ей. 
Ппаєе гесоєпійоп апа сіа5війсайоп 
0015 аге апа/угед. СІазвійсацоп ої ітаєе5 
геіаїе5 Ше прис біотеайїса! ітаге (о а рагії- 
сиїаг сіа85 ої раоїіоєїе5. Пі аддаїйоп, Фе 
апаЇузі5 ої піеап5 ої рагаПепгайоп ої Фе 
ргосе85 ої сеї писіві сопіоиг 5еіесйоп іп 
ууідевсгееп ітаєе5 из5іпе, бгарпіс ргосе55ог5 
ууаз сагтіед оці. Її 15 5помтп ас Бе рагаПе- 
П7ацйоп ої саїсиіайопя ассеіегаїез Ше рго- 
се55 ої Фаєпобі8. Ап ітаєє апаЇузі5 пподціе 
15 деуеіоред ШФаї регіогіа5 Ше ГоПом/тпе 
Гипсопз8: пріоадіпя а 5їасК ої іприї ітаєе5 
іо тетогу, рге-ргосез5іпе, ППегіпе, пріо- 
афтпя ргедейпеа райегп5, райегп 5еагспіпе 
Гог обіесів апа опгриштпя Фе па! гезціїв. 
Ап аїбогійт їог паїсго-обіесі райеги 
5еагсп 15 демеїоред уупісі соп5і5і5 ої Ше 


їоПоміпє 5(ерз: сей писіві сопіошг 5еіес- 
йоп, іфдепийсайоп ої (Фгее роїіпі5 іп сеї! 
писіеі сопіоит5, саїсціацоп ої айппіу соеї- 
Псіепі5 апа саїсиіайоп ої Фе Пієрезі 5іті- 
Іагігу ої ітаєе8. ТБе Чагабазе ої Бі5(01051- 
са! апа суіїоіовіса! іпаєе5 ої ргесапсегоце 
апа сапсегоцй5 сопдаїйоп5 ої Фе пагатагу 
5Їапд «ВРСТІ2100» уаз ц5еа Їїог Фе ехрегі- 
тепі5 ап 1( ууа5 аїзо демеїоред Бу Фе 
ашТПог8. А 50Йїмаге плодиіе м аз деуеіоред 
апа (е5іед ціп суїоіоєіса! ітаєе5 бот Ше 
ВРСІ2100 Аакабазе. Тпе демеїоред поаціе 
15 а сотропепі ої НІАМУ айіотатеа пісго- 
зсору зузіет. ТПе обіаїпед гезиіїв 5ПОму 
фас Фе еййсіепсу ої Фе рагаПепгацоп рго- 
се55 15 ароці 2090 Їог сукоїобісаї ітаєе5 
їгот Фе ВРСТ2100 4агабаєзе. 
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